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ABSTRACT 

The synchrotron radio maps of supernova remnants (SNRs) in uniform interstellar 
medium and interstellar magnetic field (ISMF) are analysed, allowing different 'sensi- 
tivity' of injection efficiency to the shock obliquity. The very-high energy 7-ray maps 
due to inverse Compton process are also synthesized. The properties of images in these 
different wavelength bands are compared, with particuh emphasis on the location of 
the bright limbs in bilateral SNRs. Recent H.E.S.S. observations of SN 1006 show 
that the radio and IC 7-ray limbs coincide, and we found that this may happen if: i) 
injection is isotropic but the variation of the maximum energy of electrons is rather 
quick to compensate for differences in magnetic field; ii) obliquity dependence of injec- 
tion (either quasi-parallel or quasi-perpendicular) and the electron maximum energy 
is strong enough to dominate magnetic field variation. In the latter case, the obliquity 
dependences of the injection and the maximum energy should not be opposite. We 
argue that the position of the limbs alone and even their coincidence in radio, X-rays 
and 7-rays, as it is discovered by H.E.S.S. in SN 1006, cannot be conclusive about the 
dependence of the electron injection efficiency, the compression/ amplification of ISMF 
and the electron maximum energy on the obliquity angle. 
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1 INTRODUCTION 

The observation of the supernova remnants (SNRs) in 
very-high energy (VHE) 7-rays by H.E.S.S. and MAGIC 
experiments is an important step toward understanding 
the nature of the Galactic cosmic rays and kinematics 
of charged particles and magnetic field in vicinity of the 
strong nonrelativistic shocks. However, the spectral analy- 
sis of multi-wavelenght data allows both for leptonic and 
hadr onic origin of VHE 7-ra y emission (e.g. RX J17 13.7- 
3946: iBerezhko fc Volkl feOQrj) . I Aharonian et all (|2007l )h In 
this context, the broad-band fitting of the spectrum of the 
nonthermal emission from SNRs is one of the hot topics in 
present studies of SNRs. At the same time, another very 
important source of scientific information, the distribution 
of the surface brightness, is not in great demand. There are 
just some discussions emphasyzing that observed correla- 
tions of brightness in radio, X-rays and 7-rays may be con- 
sidered to favor electrons to be responsible for VHE emis- 
sion in RX J1713 . 7-394 6 Vela Jr. and some other SNRs (e.g. 
lAharonian et all (|2006H . |Plagal (|2008l )V However, should the 
patterns of surface brightness in radio, X-rays and 7-rays 



realy correlate if the VHE 7-radiation originates from elec- 
trons? What should be the limitations for theory once ob- 
served patterns are really quite similar, especially in sym- 
metrical bilateral SNRs, like in SN 1006 (H.E.S.S. Source of 
the Month, August 2008). 

Another key issue for particle kinetics is the 3-D mor- 
phology of bilateral SNRs in general and SN 1006 partic- 
ularly. Is it polar-cap or barrel-like? The answer of this 
question is strongly related to the model of injection (quasi- 
parallel in the former and isotropic or quasi-perpendicular 
in the latter case) , giving therefore an important hint for ac- 
celeration theory. The properties of brightness distribution 
m ay be the most conclusi ve issue in this task (e.g. criterion 
oflRothenflug et alj (|2004l ) , azimuthal profiles comparison in 
iPetruk et al. (200$). 



An experimental investigation of SNR images have to 
be complemented with theoretical modelling of SNR maps 
in different energy domains. Radio and X-ray synchrotron 
images in the uniform interstellar medium (ISM) and the 
uniform i n terste llar magnetic field (ISMF) are modeled by 
iRevnoldsl (|l998h . The role of gradient of ISM density and 
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ISMF strength on radio m orphology of SNRs are stud- 
ied by lOrlando et all l|2007t) . These papers bases on the 
classical MHD and assumes unmodified shocks. Studies on 
nonthermal images of SNRs with non-l i near acceleration 
theory undergo development (|Lee et al.l l200Sf ). The pro- 
files of the synchrotr on brightness in such SNRs are sub- 
ject of investigation in [Ellison fc Cassam-Chena'J (|2005h and 
ICassam-Chena'f et al.l ( 200of ). 

In the present paper, we present for the first time the 
inverse-Compton 7-ray images of SNRs in uni form ISM 
and IS MF produced on the basis of the model of iRevnoldsl 
( 1998). In addition to this model, we allow for different 'sen- 
sitivity' of injection efficiency to the shock obliquity like it 
is apparent in numerical results of lEllison et al.l (|l995T ). The 
synthesized maps are compared with the radio ones. Some 
consequencies for origin of VHE emission of SNRs and elec- 
tron injection scenario are drawn. 



line of 
sight 




2 MODEL 

We consider SNR in uniform ISM and uniform ISMF. At 
the shock, the energy spectrum of electrons is taken as 
N(E) = KE~ S exp (— _E/_E max ), i5 max is the maximum en- 
ergy o f electrons , s = 2 is used throughout of this paper. We 
follow lRevnold"sl (1 19981 ) in calculation of the evolution of the 
magnet i c field and relativistic electrons (see details also in 
IPetrukl (|2006l ). |Petruk fc Beshlevl (|2008h ). The compression 
factor for ISMF era increases from unity at parallel shock to 
4 at perpendicular one. The fiducial energy at parallel shock, 
which is responsible for the ' sensitivity ' of re lativistic elec- 
trons to the radiative losses (Reynolds 1998) and which is 
used in IC images is set to -E ma x- The synchrotron losses are 
considered as the dominant channel for the radiative losses 
of relativistic electrons. We assume that K is constant in 
time; eventual evolution of K affects the radial thickness of 
rims and does not modify the mai n features of the surface 
brightness pattern (lRevnolds!ll998l ). 

Electrons emitting IC photons have energies E ~ -E max . 
Like K, E ma ,x is assumed to be constant in time. Its possible 
variation in time does not change the pattern of IC bright- 
ness and leads to effects similar to those originating from the 
time dependence of K . Namely, features in IC images have to 
be radially thicker if -E max decreases with time (i.e. increases 
with the shock velocity): since -Emax was larger at previous 
times, there are more electrons in the SNR interior able to 
emit IC photons at the present time. If i? max increases with 
time (i.e. decreases with the shock velocity) then maxima in 
brig htness are expec ted to be radially thinner. 

IRevnoldsl l|l998l) considered three models for injection: 
quasi-parallel, isotropic and quasi-perpendicular. The pat- 
tern of the radio surface brightness distribution in the 
case of the quasi-perpendicular injection is quite similar to 
the isotropic injection case, t hough with different contrasts 
l|Fulbright fc Reynolds! I1990I; lOrlando et al |[2007h . The nu- 
merical calculations of Ellison et alj 1 19951 ) show that the 
obliquity dependence of the injection efficiency <; (a fraction 
of accelereted electrons) may be either flatter or steeper than 



in the classic quasi-parallel case (? 



' O where B Q is 



the obliquity angle, the angle between the ISMF and the 
norm al to the s hock, Fig. [TJ. In order to be more general 
than IRevnoldsl (| 19981 ). we allow the injection efficiency to 



Figure 1. Geometry of the task. The obliquity angle O , the 
aspect angle O and the azimuth angle ip are shown. ISMF B is 
chosen to be parallel to the XOZ plane. 



vary with obliquity angle with different 'sensitivity' which 
is given by the parameter Ok: 



?(6o) = ?nexp(-(e /e K ) 2 ) 



(1) 



where q| is the efficiency for the parallel sh ock. This ex- 
pressi on restores approximately the results of lEllison et al.l 
(1995) with 8k = t/9 -j- 7r/4. The classic quasi-parallel in- 
jection may be approximated with Ok = 7r/6. Isotropical 
injection assumes Ok = 00, but the values Ok > 27r pro- 
duces almost the same images as Ok = 00 because the range 
for obliquity angle is < o < 7r/2. 

We consider also quasi-perpendicular injection: 

c(Oo) = ?|| exp (-((©□ - tt/2)/© k ) 2 ) • (2) 

In the most cases presented here, -E max is assumed to be 
constant over SNR surface; this c hoice allows us t o clearly 
see the role of other parameters. Reynolds (|l998l ) consid- 
ered loss-limited, time-limited and escape-limited models for 
£ ma x- In all cases, except of the loss- limited one with the 
level of turbulence comparable with the Bohm limit, i? max 
should grow with increase of ©o (|Revnoldslll9"98l ). We model 
the role of possible increase of -E max with obliquity with a 
simple parameterization 



-E ma x(0 o ) = E m&K n exp 



(-((8 



■ 7t/2)/Oe 



(3) 



where ©e is a parameter, -E max || the maximum energy at 
parallel shock. This formula, with different values of ©e, is 
able to restore approximately different cases considered by 
IRevnoldsl (|l998l) . 

The surface brightness is calculated integrating emis- 
sivities along the line of sight within SNR. The synchrotron 
emissivity at some radio frequency is q syc h oc KB^ 3+1 '' 2 , 
B is the strength of magnetic field. The 7-ray emissivity of 
electrons due to inverse Compton process is calculated as 



N(E)p ic (E,e)dE 



(4) 
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Figure 2. Radio images of SNR for an aspect angle 4>o = 90° and different 0k: 7r/12 (a), 7r/6 (b), 7r/4 (c), tt/2 (d), ir (e), 2tt (f). 
Ambient magnetic field is oriented along the horizontal axis. Hereafter, the increment in brightness is AS = O.lSmax- 



where e is the photon energy. The spectral distribution p lc 
of radiation power of a "single" electron in a black-body 
photon field with temperature T is 



Pic(7,e) 



77 2ic(r? c ,?7o) 



where 7 is Lorenz factor of electron, e c = kT, 



Vc = 



(m c c 2 ) 2 ' 7,0 4~fin c c 2 (jm c c 2 - e) ' 



(•>) 



(6) 



m e , e, c, 7i, k hav e their typica l meaning. 2i c (?7 c ,?7 ) may be 
approximated as (|Petrukll2008h 



l ic (Vc,Vo) » y?7c exp 



+2r/ exp 
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3r? c 



(7) 



This approximation is quite accurate, it represents Ii c in 
any regime, from Thomson to extreme Klein-Nishina. The 
maximum of s pectral distrib ution Pi c (e) for electrons with 
energy E is at (|Petrukll2008l ) 



r c = 



ie c E 
(m e c 2 ) 2 



(8) 



All IC images in the present paper (except of that on Fig. llOf) 
are calculated for the initial photon field with T = 2.75 K 
and for the 7-ray photon energy e = 0.1e ma x(i?max) that is 
for example e = 0.3 TeV for £ max = 30TeV. 



3 RESULTS 

3.1 Synchrotron radio images 

We stress that all figures in the present paper have been 
computed using complete MHD model. 

Let us define an aspect angle (j> as an angle between 
interstellar magnetic field and the line of sight (Fig. [TJ. It 
is shown that the azimuthal variation of the radio surface 
brightness S g at a given radius of projection g, in SNR which 
is not centrally brightened, is mostly determined by the vari- 
ations of the magnetic field compression (and /or amplifica- 
tion) (jb and the electron injection efficiency <; (|Petruk et al.l 
120091) : 

S e (<p) x ?(6 D> sff(¥>,0a)) o-B(O o , c ff(v5,0 o )) <S+1)/2 (9) 

where tp is the azimuthal angle. The effective obliquity angle 
Oo.eff is related to ip and <j> as 



cos 9o, c ff (if, <f>o) = cos if sin O 



(10) 



here, the azimuth angle 93 is measured from the direction of 
ISMF in the plane of the sky (Fig.[T]). 

Fig. [2] shows how 0k affects a radio image of SNR. 
Complete MHD simulations are in agreement with the ap- 
proximate formula ([9]). First, we note that smooth increase 
of Ok results in transition from the 3-D polar-cap model of 
SNR to the 3-D barrel-like one. This is also visible on Fig. [3] 
where ISMF is directed toward observer. Namely, increase 
of Ok change the visual morphology from centrally-bright 
to shell-like. 

There are three names for a class of SNRs which 



l|Kesteven & Caswelll 


1987), 'bipolar' 


199d) and 'bilateral' 


Gacnslcr 1998) 



on the base of 2-D visual morphology. It is interesting that 
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Figure 3. Profiles of the radio surface brightness for an aspect 
angle <f> = 0° (the radial profile of brightness is the same for any 
azimuth). B K is tt/12 (line 1), tt/6 (line 2), tt/4 (line 3), tt/2 (line 
4), 7T (line 5), 2n (line 6). 

the first two names reflects de facto the two different con- 
ceptions of SNRs in 3-D. 

Fig. [5] also shows that an assumption about orienta- 
tion of ISMF leads to limitation of possible injection model. 
Ambient magnetic field in all images on Fig.[2]is along hori- 
zontal axis. Thus, if one consider the polar-cap scenario for 
bilateral SNR (ISMF is along axis which crosses two limbs) 
then one should consider the injection model which strongly 
depends on the obliquity (8k < tt/6, Fig. Oi,b). Instead, if 
the barrel is the preferable model (ISMF is parallel to the 
symmetry axis between two limbs) then the injection effi- 
ciency should be almost independent of obliquity (O > n, 
Fiff. j2fe,f), or prefer quasiperpendicular shocks. 

iGaenslerl (|l998l ) measured the angle tp between the sym- 
metry axis in 17 'clearly' bilateral SNRs and the Galactic 
plane. Axes are more or less aligned with the Galactic plane 
in 12 SNRs (ip < 30°), 2 SNRs have tp « 45° and 3 SNRs is 
almost perpendicular {ip > 60°). If we assume that ISMF is 
parallel to the plane of Galaxy then most of bilateral SNRs 
should be 3-D barrels preffering thus isotropic (or quasiper- 
pendicular) injection. 

An interesting feature appears on images for Ok = 
7r/4-f7r/2 (Fig.[2p,d). Namely, SNR has 'quadrilateral' mor- 
phology. With increasing of obliquity, the injection efficiency 
decreases while the compression factor of ISMF icreases. The 
variation of injection ?(©o) dominates o"b(©o) for Ok < n/6. 
If 6k > 7T (injection is almost isotropic) then <tb(©o) plays 
the main role in azimuthal variation of the radio surface 
brightness. In the intermediate range of Ok, the signifi- 
cance of the two variations are comparable leading there- 
fore to azimuthal migration of the brightness maxima in the 
modelled images. There is no 'quadrilateral' SNR reported 
in the literature. If there is no such SNR at all, the range 
Ok — 7r/4-^7r/2 may be excluded. However, we stress that a 
complete statistical study of the morphology of radio SNRs 
would be needed to definitly asses the lack of quadrilateral 
SNRifl 

The visual morphology of SNR is different for different 
aspect angles. Fig. [3] shows SNR images for quasi-parallel 
injection with Ok = ir/12 (upper panel) and for isotropic 
injection (Ok = 27r, lower panel). We may expect that 

1 G338. 3-0.0 could be an example of quadrilateral SNR 



ISMF may have different orientation versus observer in var- 
ious SNRs. If quasi-parallel injection is not a rare exception 
then the polar-cap SNRs should be projected in a different 
way and we may expect to observe not only 'bipolar' SNRs 
(Fig. |4p,d) but also SNRs wit h one or two radio eyes within 
thermal X-ray rim (Fig. lib.b). iFulbright fc Reynolds! (|l990h 
developed statistically this thought and showed that the 
quasi-parallel injection model would be unlikely, but again, 
we would need a comple te study to verify this statem enfl 
Statistical arguments of iFulbright fc Reynolds! {l990) may 
be affected by the fact that centrally-bright radio SNRs 
(lines 1-2 on Fig. [3J) are expected to be fainter than bilateral 
or circular SNRs with the same characteristics (lines 4-6 on 
Fig. [3}: it could be that most of the centrally-peaked SNRs 
may not be observable. 



3.2 IC 7-ray images 

Let us consider first the case when the maximum energy 
of electrons is constant over SNR surface; this allows us to 
clearly see the role of the injection efficiency and magentic 
field variations. 

Synthesized IC 7-ray images of SNRs are presented on 
Fig. for different aspect angles. These images assume al- 
most isotropic injection (Ok = 2n) and should be compared 
with radio maps on the lower panel of Fig. [3] The compo- 
nent of ISMF which is perpendicular to the line of sight is 
along horizontal axis on all images. An important difference 
is prominent from these two figures. Namely, the two lobes 
develop with increasing of <f> in both radio and 7-rays. How- 
ever, their location in respect to ISMF is opposite. The line 
conecting two maxima in radio is perpendicular to ISMF 
while it is parallel to ISMF on IC images (cf. Fig [SJi and 
Fig Si). 

The reason of this effect is the following. For assumed 
isotropic injection, the azimuthal variation of the radio 
brightness is determined only by the dependence <tb on 
obliquity (the azimuth angle equals to the obliquity angle 
for <j) = tt/2). Electrons emitting VHE 7-rays have energies 
E ~ -Emax and experience substantial radiative losses (this 
effect is negligible for radio emitting electrons). Magnetic 
field does not appear directly in the formulae for IC emis- 
sion, but it affects the downstream distribution of relativistic 
electrons emitting IC 7-rays. The larger post-shock magnetic 
field the larger radiative losses. The downstream distribution 
of IC-emitting electrons is therefore steeper where magnetic 
field is stronger. This leads to lower IC brightness in SNR 
regions with larger magnetic field (while radio brightness 
increases there because of proportionality to B 3 ^ 2 ). 

In VHE 7-ray image of SN 1006 recently reported by 
H.E.S.S. collaboration (H.E.S.S. Source of the Month, Au- 
gust 2008), the two maxima coincide in location with limbs 
in radio and nonthermal X-rays. This fact, in view of the 
'limb-inverse' property, could be considered as argument 
against the leptonic origin of 7-ray emission in SN 1006 
(if injection is isotropic). However, these IC images are ob- 
tained under assumption that E max does not vary over SNR 



2 G311. 5-0.3 and G337.2-0.7 could be examples of SNRs with two 
radio 'eyes' 
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Figure 4. Radio images of SNR for different aspect angles <f> : 0° (a,e), 30° (b,f), 60° (c,g), 90° (d,h). 6 K = 7r/12 (upper panel), 6 K = 2vr 
(lower panel). Component of the ambient magnetic field which is perpendicular to the line of sight, is oriented along the horizontal axis. 
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Figure 5. IC 7-ray images of SNR. Isotropic injection, i? ma x is constant over SNR surface. Aspect angles <j> : 0° (a), 30° (b), 60° (c), 
90° (d). Component of the ambient magnetic field which is perpendicular to the line of sight, is oriented along the horizontal axis. 



surface. If i? max is high enough at regions with large mag- 
netic field (at perpendicular shock), then the 'limb-inverse' 
effect may be less prominent or even might not be important 
(see below). 

In case if injection strongly prefers parallel shocks (limbs 
in SN 1006 are polar caps), the dependence c(©o) might 
dominate ctb(©o)- The maxima of brightness in radio and IC 
7-rays are therefore located at the same regions of SNR pro- 
jection (Fig.[6j to be compared with Fig.|4^,,d), in agreement 
with the Chandra and H.E.S.S. observations of SN 1006. 

The role of intermediate values Qk for injection which 
prefers parallel shock, Eq. ([TJ, on profiles of IC brightness 
is shown on Fig. [7] Increase of the sencitivity of injection 
to the obliquity leads to radially thinner and more contrast 
features. 

If injection prefers perpendicular shock, Eq. ([2]), its in- 
crease in the regions of larger magnetic field may compen- 
sate the lack of 7-ray emitting electrons. In that case, the 
position of limbs coincide in radio and IC 7-rays if the de- 



pendence ?(© ) is strong enough (Fig.|8p,d). In the range of 
intermediate 0k, the quadrilateral morphology appears also 
in models of IC 7-rays (Fig.|8t), as an intermediate morphol- 
ogy between those on Fig. [5jl and Fig. |HJl. (The contrast of 
maxima in the image of quadrilateral SNR is so small that 
this feature may probably not be observable.) 

Note that the quasi-perpendicular injection model leads 
to radio images similar to those in the isotropic injectio n 
case, cf. Fig.|5K,b and Fig.[2f (see also lOrlando et all l|2007h ). 
because magnetic field and injection efficiency increase at 
perpendicular shocks both leading to larger synchrotron 
emission. In contrast, there is a lack of IC radiating elec- 
trons around perpendicular shocks which may or may not 
(depending on 6k in ([2])) compensate it. Thus IC images 
involving the quasi-perpendicular injection may radically 
differ from those with isotropic injection, cf. Fig. |8jl and 
Fig.EJi. 

The obliquity variation of the electron maximum energy 
is an additional factor affecting the IC 7-ray brightness in 
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Figure 6. IC 7-ray images of SNR. Quasi-parallel injection JIJ 
with © K = 7r /4, -E m ax(©o) = const. Aspect angles <j> : 0° (a), 90° 
(b). In the latter, ISMF is along the horizontal axis. 



SNRs. Actually, iRothenflug et all (|2004 ) have shown that 
the cut-off frequency increases at radio limbs of SN 1006 
that may (partially) be due to larger -E max there. Therefore 
£ mal is expected to be largest in this SNR at the perpen- 
dicular shock (at equatorial belt) if injection is isotropic or 
quasi-perpendicular or at the parallel shock (at polar caps) if 
in jection is quasi-p arallel. In the latter case, the calculations 
of I Reynolds! (|l998T ) suggest that the only possible model for 
£max in SN 1006 should be loss-limited one in the Bohm 
limit. 

The role of -Emax increasing with obliquity, Eq. Q, is 
shown on Fig. [9] The 'limb-inverse' property may not be 
important and the limbs may coincide in radio, X-rays and 
IC 7-rays even for the isotropic injection if the maximum 
energy is large enough at perpendicular shocks to provide 
energetical electrons in despite of radiative losses (Fig. [9p, 
cf. with Fig. UJi and Fig. [SJl). Note also that the limbs are 
thicker in this case, because of the more effective radiative 
losses at perpendicular shock (due to larger ISMF compres- 
sion), comparing to limbs if they are at parallel shock. 

The dependence of E max on O may also cause splitting 
and rotation of IC limbs in case of the quasi-parallel injection 
(Fig. |5Ji, cf. with Fig. [Up) or the quasi-perpendicular one. 
There is a possibility for quadrilateral SNRs to appear in 
7-rays due to the interplay between dependences E m ax(0 o ), 
?(0 O ) and cr B (0 o ) (Fig.Ek,d). 

All above IC images are calculated for the photon en- 
ergy e = 0.1e ma x(£'max)- The pattern of the 7-ray surface 
brightness remain almost the same with increasing of the 
photon energy, though regions of maximum brightness be- 
come radially thinner and also contrasts change (Fig. I10[) . 
This is because electrons which contribute most of emission 
at larger photon energy experience higher radiative losses 
and therefore the downstream distribution of these electrons 
are steeper. 

To the end, the main properties of IC surface brightness 
may simply be derived from the approximate analytical for- 
mula for the azimuthal variation of IC surface brightness 
S e {ip\4> ,e) of the adiabatic SNR in uniform ISM and uni- 
form ISMF (Appendix): 



S g (ip) oc ?(e o , e ff)exp - 



F n- 1 - 5CT B(e o , cff ) 2 B m /2B fi 



J raax, || 



(11) 




where E m oc e 1 ^ 2 , Eq. (|A8[) , g = g/R < 1, g is the distance 
from the center of SNR projection. This formula may not be 



Figure 7. Profiles of the IC surface brightness along X-axis for 
the aspect angle <p = 0° (the radial profile of brightness is the 
same for any azimuth; to be compared with Fig. [3j and <f>o = 90° 
(ISMF is along the horizontal axis). Dependence of injection is 
given by JTJ with 0k (from below): 7r/12, 7r/6, it/4, tt/2, tt, 2tt, 
00. -Emax is constant over SNR surface. 



used for SNR which is centrally-bright in 7-rays and is valid 
for g/R larger than ~ 0.9. 



4 CONCLUSIONS 

In the present paper, we analyse the synchrotron radio and 
the inverse-Compton 7-ray images of Sedov SNRs synthe- 
sized on the base of the lRevnoldsi ([1998) model. IeU ison et al.l 
(1995) have shown that the dependence of efficiency of in- 
jection c on obliquity angle O may differ from commonly 
used expression in quasi-parellel case. We therefore param- 
eterise the dependence ?(0 O ) as it is given by Eq. JT]). It 
is shown that the variation of the parameter 0k provide 
smooth transition from polar-cap (0k < t/6) to barrel-like 
(0k > tt) models of SNR and that assumed orientation of 
ISMF should be related to a certain injection model. Some 
constraints on injection models which follow from morpho- 
logical considerations are pointed out. The azimuthal varia- 
tion of radio brightness is mostly due to variations of ? and 
(Jb, in agreement with the approximate formula ©. 

Theoretical 7-ray images of SNR due to the inverse 
Compton effect are reported for the first time. We analyse 
properties of these images and compare them with corre- 
sponding radio maps of SNRs. The azimuthal variation of 
IC brightness is mostly determined by variations of era 
and -E max , in agreement with the approximate formula (|ll|l 
derived in the Appendix. 

In case if E max is constant over the SNR surface, we 
found an opposite behaviour of azimuthal variation of sur- 
face brightness in radio and IC 7-rays, in case if injection is 
isotropic and the aspect angle is larger than ~ 60°. Namely, 
the line crossing the two limbs in radio are perpendicular to 
the ISMF while they are parallel in IC 7-rays. In particu- 
lar, bright radio limbs correspond to dark IC areas, in dis- 
agreement with X-ray and H.E.S.S. observations of SN 1006. 
This happens because IC image is affected by large radia- 
tive losses of emitting electrons behind perpendicular shock 
while the larger magnetic field increases the radio brightness 
there. Variation of E max over SNR surface may (to some ex- 
tent) hide this effect. The maximum energy should increase 
with obliquity in this case. 
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Figure 8. Radio (a,b) and IC 7-ray images (c,d) of SNR for 
<j>o = 90°. Quasi-perpendicular injection {2} with 0k = t/2 (a,c) 
and ©k = t/6 (b,d) (to be compared with Fig.|4jl and Fig. |5jl). 
Bmax is constant over SNR surface. 



In case of the polar-cap model of SNR (quasi-parallel 
injection), the maxima in surface brightness are expected to 
coincide in radio and IC 7-rays (in agreement with H.E.S.S. 
observation of SN 1006), unless increase of -E max with obliq- 
uity will be very strong, which is unlikely in case of SN 1006 
because the cut-off frequency is larger at limbs which are at 
parallel shock in this injection model. 

Limbs may also coincide in case of the quasi- 
perpendicular injection, if the lack of electrons (due to ra- 
diative losses) in the regions of large magnetic field is com- 
pensated by the strong enough increase of ? and/or _E max 
with o . 

Isotropic compression/amplification of ISMF on the 
shock (i.e. independent of the shock obliquity), like it could 
be under highly effective acceleration, may also be responsi- 
ble for the same position of limbs in radio and in IC 7-rays, 
for the quasi-parallel or quasi-perpendicular injection sce- 
narios. In this case the dependence of -E max (0 o ) have to 
follow variation ?(0 O ), namely, to be largest (smallest) at 
parallel shock for quasi-parallel (quasi-perpendicular) injec- 
tion, otherwise the morphology of SNR in IC 7-rays may 
differ from the radio one. 

We conclude that the location the 7-ray limbs ver- 
sus radio and X-ray ones, recently discovered by H.E.S.S. 
in SN 1006, cannot be conclusive about the actual de- 
pendence of the electron injection efficiency, the compres- 
sion/amplification of ISMF and the electron maximum en- 
ergy on the obliquity angle in this SNR. Detailed features of 
the SNR maps in different wavebands should be considered 
for this purpose. 

The interplay between dependences ?(Oo), ob(Oo) and 
(0 O ) may cause the quadrilateral morphology in SNR 
models, due to splitting of maxima in surface brightness. 
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Figure 9. IC 7-ray images of SNR for if>o = 90° and E max 
increasing with obliquity, Eq. Q with ©e = 7r/2 (a,c) and 
©E = 7r /4 (b,d). Isotropic injection (a,b), to be compared with 
Fig.[5jl; quasi-parallel injection with ©k = 7r/4 (c,d), to be com- 
pared with Fig. [6b. 
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Figure 10. The same as Fig. Q5J1, for 10 times larger photon 
energy, e = Emax(Emai). 

Absence of quadrilateral SNRs in IC 7-rays, if revealed ob- 
servationally, may results in limitations on 0k and 0e. 

The detailed characterictics of features on IC image 
(e.g. thickness of rim) depend on the photon energy. They 
are radially thinner at larger photon energies, as expected. 



ACKNOWLEDGMENTS 

OP acknowledge Osservatorio Astronomico di Palermo for 
hospitality. The work of OP was partially supported by the 
program 'Kosmomikrofizyka' of National Academy of Sci- 
ences (Ukraine). FB, SO and OP acknowledge Consorzio 
COMETA under the PI2S2 Project, a project co-funded 
by the Italian Ministry of University and Research (MIUR) 



8 Petruk O. et al. 



within the Piano Operativo Nazionale "Ricerca Scientifica, 
Sviluppo Tecnologico, Alta Formazione' (PON 2000-2006). 



REFERENCES 

Aharonian, F. et al., 2006, A&A 449, 223 
Aharonian, F. et al., 2007, A&A 464, 235 
Berezhko, E. G., & Volk, H. J. 2006, A&A 451, 981 
Cassam-Chenai G., Decourchelle A., Ballet J., Ellison D. 

C, 2005, A&A 443, 955 
Ellison D. C, Baring M. G., Jones F. C, 1995 ApJ, 453, 

873 

Ellison D. & Cassam-Chenai G., 2005, ApJ, 632, 920 
Gaensler B. M., 1998, ApJ, 493, 781 
Fulbright M. S., & Reynolds S. P., 1990, ApJ, 357, 591, 
Kesteven M. J. & Caswell J. L., 1987, A&A, 183, 118 
Lee S.-H., Kamae T., Ellison D. C, 2008, ApJ, 686, 325 
Orlando S., Bocchino F., Reale F., Peres G., & Petruk O., 

2007, A&A, 470, 927, 
Petruk O., 2006 A&A, 460, 375 
Petruk O., 2008, astro-ph/0807.1969 
Petruk O., Beshley V., 2008, KPCB, 
Petruk O., Dubner G., Castelletti 

Kirsch M., Miceli M., Orlando S. 

MNRAS, accepted 
Plaga R., 2008, New Astronomy, 13, 
Reynolds S. P., 1998, ApJ, 493, 375 
Rothenflug R., Ballet J., Dubner 

courchelle A., & Ferrando P., 2004, 
Schlickeiser R. Cosmic Ray Astrophysics (Springer, 2002) 



24, 159 

G., Iakubovskyi D., 
Telezhinsky I., 2009, 



73 

G., Giacani E.. 
A&A, 425, 121 



De- 



APPENDIX A: APPROXIMATE ANALYTICAL 
FORMULA FOR THE AZIMUTHAL 
VARIATION OF THE IC 7- RAY SURFACE 
BRIGHTNESS IN SEDOV SNR 

An approximate formula for azimuthal variation of the IC 7- 
ray surface brightness allows one to avoid detailed numerical 
simulations and may be useful if approximate estimation for 
the variation is reasonable. It gives deeper insight in the 
main factors determining the azimuthal behavior of the IC 
surface brightness in SNRs. 

Let the energy of relativistic electrons is E in a given 
fluid element at present time. Their energy was Ei at the 
time this element was shocked. These two energies are re- 
lated as 



E = EiEadSn 



(Al) 



where £ a( j accounts for the adiabatic losses and £ ra d for the 
radiati ve losses. There are app roximations valid close to the 
shock l|Petruk fc Beshlevll2008l ): 

£ ad ~ a, £ rad ^ a 5CT B E / 2E f,ll (A2) 

where a — a/R, a is Lagrangian coordinate of the fluid el- 
ement, i?fj is the fiducial energy for parallel shock. The 
downstream evolution of K in a Sedov SNR is 



K oc q{Q )K(a). 
With the approximations 



be written in the model of iRevnoldsl (|l998l ) as 



(A3) 

the distrib ution N(E) may 



N(E, e o ) oc s(0 o )K(a)E- s exp 



En -*(E,e o) 



E„ 



11^(60) 



where 



ip(E,& ) = 1 + 



5o- B (e ) 2 g 

2E t 11 



(A4) 



(A5) 



and the obliquity variation of the maximum energy of elec- 
trons is given by £ m „ = Smax.H-^Qo)- 

Electrons with Lorentz factor 7 emit most of their IC 
radiation in photons wit h energy e m . Let us use the 'delta- 
function approximation ' l|Petrukll2008h : 



Pic (7, 



: Pm(7)5(e- 



), Pm(7) = Pic( r r,s)de. (A6) 







In the Thomson li mit, which is valid for SNRs in most cases, 
emf-y) « 4fc?V (IPetrukl 120081 ) and p m [i) = (4/3)c<j T w7 2 
( Schlickeiser 2Q0l), T and lj are the temperature and the en- 
ergy density of initial black-body photons, <7t is the Thom- 
son cross-section. 

Substitution @ with (|A6|) yields 



co"Twm e c 2 e 1 '' 2 



12e; 



3/2 



N(E m ) 



where 



E m = 



m c cV /2 
2(feT)V2 



(A7) 



(A8) 



is the energy of electrons which give maximum contribution 
to IC emission at photons with energy e. 

Let us consider the azimuthal profile of the IC 7-ray 
brightness S e at a given radius g from the centre of the 
SNR projection. 

The obliquity angle O is different for each radial sec- 
tor of 3-D object. It is determined, for any position within 
SNR, by the set (ip,f / 'q,4> ). Integration along the line of 
sight gathers information from different radial sectors, with 
different obliquities. Let us determine the 'effective' obliq- 
uity angle by the relation 



6 ,eff(¥>, <t>o) = ©o(<p, 1, 4>o)- 



(A9) 



Actually, O , ff for a given azimuth equals to the obliquity 
angle for a sector with the same azimuth lying in the plane 
of the sky (i.e. in the plane being perpendicular to the line 
of sight and containing the center of SNR). o varies around 
o ,eff during integration along the line of sight. The closer 
g to the edge of SNR projection the smaller the range for 
variation of O and more accurate is our approximation. 

The surface brightness of SNR projection at distance g 
from the center and at azimuth if is 



S{q, V )=2\ q ic (a) _ 
8(e) V r 



rrada 



(A10) 



where f a is the derivative of r(a) in respect to a. The az- 
imuthal variation of the IC brightness for fixed g is approx- 
imately 
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S e oc c(9o,eff) exp 




) 



(All) 



x 



a(e) \/r 2 - Q 2 



rr^da 



exp 




) 



If q — *• 1 then cl(q) — > 1. Thus, the exponent in the integral 
is roughly unity because a(g) < a < 1 and a(g) < g < 1. 
The integral in (| Al 1 1) is therefore roughly the same for any 
azimuthal angle ip. The azimuthal variation of the IC 7-ray 
brightness S e (ip; (f> , e) is thus determined mostly just by 



with E m given by Eq. ()A8|) . i.e. S e depends in this approx- 
imation on the temperature T of the seed black-body pho- 
tons and the energy e of observed 7-photons. The relation 
between the azimuthal angle if, the obliquity angle o ,eff 
and the aspect angle <f> is as simple as 

cos O 0j efi (<p, </> ) = cos ip sin 4>o (A13) 

for the azimuth angle <p measured from the direction of ISMF 
in the plane of the sky. 

The approximation (|A12() may be used for g larger than 
~ 0.9. Like Eq. ©, Eq. (|A12|I does not give correct profiles 
in the case of centrally-bright SNRs, i.e. when Ok < 7r/4 in 
© and <j!> < 30°. 




